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Abstract 

We derive the two-loop Bethe ansatz for the sl(2) twist operator sector of N = 4 gauge 
theory directly from the field theory. We then analyze a recently proposed perturbative 
asymptotic all-loop Bethe ansatz in the limit of large spacetime spin at large but finite 
twist, and find a novel all-loop scaling function. This function obeys the Kotikov-Lipatov 
transcendentality principle and does not depend on the twist. Under the assumption that 
one may extrapolate back to leading twist, our result yields an all-loop prediction for the 
large-spin anomalous dimensions of twist-two operators. The latter also appears as an 
undetermined function in a recent conjecture of Bern, Dixon and Smirnov for the all-loop 
structure of the maximally helicity violating (MHV) n-point gluon amplitudes of M = 4 
gauge theory. This potentially establishes a direct link between the worldsheet and the 
spacetime S-matrix approach. A further assumption for the validity of our prediction is 
that perturbative BMN (Berenstein-Maldacena-Nastase) scaling does not break down at 
four loops, or beyond. We also discuss how the result gets modified if BMN scaling does 
break down. Finally, we show that our result qualitatively agrees at strong coupling with 
a prediction of string theory. 



1 Introduction and Main Results 



There is mounting evidence that planar M = 4 gauge theory might be " exactly solv- 
able" . For example, it was recently proposed that higher-loop maximally helicity vi- 
olating (MHV) n-point gluon amplitudes should be iteratively expressible through the 
(regulated) one- loop amplitudes PQ, [2]. Based on sophisticated two- loop 0,[T] and three- 
loop [2] computations, a conjecture for the all-loop MHV n-point gluon amplitudes M. n 
in 4 — 2 e dimensions was formulated in [2] : 



M. n = exp 



(1) 



Here En (e) vanishes as e — > 0, are finite constants, and M n (£e) is the (£e)- 
regulated one-loop n-point amplitude. At e = we have lima e = g 2 where g 2 is defined 



as 



^=^T = ^ ( 2 ) 



9ymN _ 
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and A is the 't Hooft coupling. Finally, the /^(e) are generated in the e — ► limit by 
the function 

oo 

f(g)=AY,9 2e f (e) (0). (3) 

This function is in fact related to the large spin anomalous dimension of so-called leading- 
twist operators in the gauge theory [3]. Alternative names are "soft" anomalous dimen- 
sion, or "cusp" anomalous dimension. The simplest representatives of Af = 4 twist 
operators are found in the sl(2) sector: 

Tr (V s Z L ) + ... . (4) 

Here T> is a light-cone covariant derivative, Z is one of the three complex scalars of the 
M = 4 model, s measures the spacetime spin, and the twist L is to equal to, in this sector, 
one of the so(6) R-charges. Leading twist is L = 2. The dots in (jlj) indicate that the true 
quantum operators are complicated mixtures of states, where the s covariant derivatives 
may act in all possible ways on the L fields Z. Mixing with multi-trace operators is 
suppressed in the planar theory. The function f(g) is obtained by considering the large 
spin s — > oo limit of the anomalous scaling dimension of the quantum operators (j3J), 
which is expected to scale logarithmically as 

A = s + /(<?) log(s) + O(s ) . (5) 

We will call f(g) of the scaling function. Note that the scaling structure in © is 

a highly non-trivial structural property of the exact finite s expression for A = A(s,g). 
Individual Feynman diagrams contributing at intermediate stages of the perturbative 
calculation of A certainly contain higher (k > 1) powers log fc (s). The one-loop 0(g 2 ) 
contribution to A was first computed in [51 [6] for all s, and indeed behaves as in (J5]). 
The 0(g A ) two-loop answer was found in [HE], an d the 0(g 6 ) three-loop one, inspired 
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by a full-fledged computation in QCD [9], in [10]. Again, the result indeed scales as in 
([5]), and the state of the art up to now has been [10] 

f{g)=±g 2 -\* 2 g* + l ^g* + .... (6) 

Fascinatingly, this agrees via ([31) with the three-loop n = 4 calculation of ([1]) by Bern, 
Dixon and Smirnov [2]. 

There is also mounting evidence that planar Af = 4 theory might be "exactly inte- 
grable". This means that the spectral problem, i.e. the spectrum of all possible scaling 
dimensions {A} of the Af = 4 gauge theory, is encoded in a Bethe ansatz. This was estab- 
lished for the complete set of possible operators of Af = 4 theory at one loop [HI [12]. It 
was then conjectured, based on two- and three-loop computations in the model's su(2) 
sector, that integrability extends to all orders in perturbation theory, and, hopefully, 
to the non-perturbative level [13]. This was subsequently backed up in various studies 
[141 [T5l \1Q\ fm [T8| [T9] , and culminated in a proposal for the asymptotic all-loop Bethe 
equations of the theory [20], further supported by [21]. Here "asymptotic" means that 
the ansatz of [20] is only expected to correctly yield the anomalous dimension up to, 
roughly, 0(g 2L ~ 2 ), where L is the number of constituent fields (not counting covariant 
derivatives) of the considered quantum operators^. Very recently a proposal was made 
for circumventing the asymptoticity restriction for the su(2) sector by relating the di- 
latation operator, whose eigenvalues are the dimensions A, to a Hubbard Hamiltonian 
[22]. It would be exciting to also "hubbardize" the sl(2) sector of twist operators. 

Let us stress that "solvability" and "integrability" are distinct concepts. The latter 
is a rather precise but narrow concept referring to the spectral problem of the gauge 
theory. It means that the spectrum is described by one-dimensional factorized scattering 
of a set of appropriate elementary excitations. Equivalently, it means that there is a 
Bethe ansatz. Exactly solving the Bethe ansatz, in a given situation, is rarely easy, and 
actually generically impossible. Integrability allows one to prove that there will never be 
a "plug-in" formula for the spectrum of Af = 4 gauge theory! 

On the other hand, "solvability" is a significantly less precise concept. Nevertheless, 
there is much evidence that in Af = 4 gauge theory many quantities beyond the scaling 
dimensions allow for a precise mathematical description. We have begun our discus- 
sion with the conjecture (CQ), which clearly contains more than just spectral information. 
Another example are coordinate space correlation functions of more than two local com- 
posite operators. Certain intriguing iterative structures were e.g. noticed in four-point 
functions some time ago [23j EH [25] . At the time of writing, the precise relation between 
the observed solvable structures and the integrable structures is somewhat reminiscent 
of the well-known paradox of the chicken and the egg. 

Recall that a three-loop st(2) Bethe ansatz for gauge theory was conjectured in [18] 
by taking inspiration from the integrable structures appearing in string theory |26j , and 
indeed reproduced ©• We will further back up the conjecture of [IE] by calculating 
in section [2] the two-loop S-matrix of the sl(2) sector directly from the field theory. 

1 In the case of twist L = 2 operators the asymptotic ansatz actually works to 0(g 6 ) instead of the 
naively expected 0(g 2 ) because of superconformal invariance [18l 120] . 
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An alternative two-loop derivation, using algebraic methods, was recently presented by 
Zwiebel [19J (actually, for the bigger sector su(l, 1|2)). 

It is amusing to note that §§§ is thus reproduced by three completely independent 
procedures ([T0].,[T%].[!Z]). none of them completely rigorous. However, the various ap- 
proaches, including their assumptions, seem to be completely independent. So ([6]) is very 
likely to be correct! 

In this paper we will apply the asymptotic all-loop Bethe ansatz of [20] in order to 
compute all further perturbative corrections to the expression Strictly speaking, the 
asymptotic ansatz does not apply to leading twist L = 2, see above. We will however 
argue that the scaling function ([5]) is universal in that it describes the behavior of the 
lowest state of any st(2) operator as long as I C s. For a very recent discussion of 
this point, on the one-loop level, see [27]. Our argument for the validity of our scaling 
function, as concerns the leading twist operators, is therefore based on two assumptions: 
(1) That it is indeed correct to pick L sufficiently large to stay in the "asymptotic" 
regime of the Bethe ansatz, while keeping L C s, and (2) that the Bethe ansatz of [20] 
indeed describes the gauge theory, for sufficiently "long" operators, at four loops and 
beyond. The first assumption is very likely to be true, while the validity of the second is, 
at the time of writing, much less clear. However, our computation might actually help 
to decide this issue, see below. 

As a highly non-trivial check of our procedure, we will prove that the anomalous 
dimension A is indeed of the expected scaling form ([5]) to all order in perturbation 
theory. We will find the scaling function to be given by the integral representation 



POO 

f(g) = Ag 2 - 16 s 4 / dta(t) 
Jo 



Ji(V2gt) 
V2gt ' 



(7) 



where the fluctuation density a(t) is determined by the solution of the integral equation 
t 



a(t) 



e* - 1 



Jl ^ 9t) -2g 2 rdt'K(V2gt,V2gt>) &(t>) 

V2gt Jo 



with the non-singular kernel 

£( t ff) = Jl(t)J Q (t') ~ Jo(t) Mf) (g) 

The functions Jo(t), Jxif) i n the above equations are standard Bessel functions. 

We have been unable to find an explicit solution of the integral equation. It would 
be quite interesting if this could be achieved. It is however straighforward to obtain 
the weak-coupling expansion of the fluctuation density a(t) by iterating (jSJ). Using ((7j) 
we then obtain the perturbative solution of the scaling function f(g). To e.g. four- loop 
order one has 

f(g) = Ag 2 -4C(2)s 4 + (4((2) 2 + 12C(4))^ 6 (10) 
- (4 C(2) 3 + 24 C(2)C(4) - 4 C(3) 2 + 50 C(6)) / + ..., 
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Using the fact that ^-functions of even argument may be expressed as products of rational 
numbers and powers of n, this may be simplified to 

/(<?) = W - I «* 9* + % 9 6 ~ - 6 - 4 C (3) 2 ) / + .... (11) 

As a further non-trivial check of our procedure, and thus the validity of (J7|), we 
shall find that f(g) obeys the Kotikov-Lipatov principle of maximal transcendent ality 
[7j, which was actually used in [TU] in order to extract, even at finite spin s, the J\f = 4 
dimensions from the QCD calculation of [9J. When applied to the large s limit, the 
principle holds that the sum over all the arguments of the products of zeta functions 
appearing as additive terms at a given loop order £ always adds up to 2£ — 2, see f fTUj) . 
and ([51]) below. 

It would be exciting if the four-loop prediction ffTTj) could be tested by a field theoretic 
computation, maybe by way of extending the results of [2] to higher order. Incidentally, 
even if field theory fails to reproduce (ITT]) , we will gain crucial knowledge on the integrable 
structure of the gauge theory, see section 13.41 The reason is that we are able to predict 
how transcendentality will break down if it breaks down. The e.g. four-loop term in the 
scaling function f(g) in (fiT]) would then get modified to 

- (^ 6 -4C(3) 2 + 8/3C(3)) 9*, (12) 

where (3 is an a priori unknown numbeJE Furthermore this number would then show 
up in the four-loop anomalous dimensions of all operators of the M = 4 theory, see 13.41 
In particular, it would manifest itself in the four-loop dimensions of operators with a 
large R-charge J, and would in fact induce a perturbative breakdown of BMN-scaling 
[28] . The argument may also be turned around: Proving that (flT!) holds as stated would 
establish that (3 = 0, and would therefore be indirect proof that BMN-scaling holds up 
to the four-loop level. 

Finally, there is a prediction from string theory [291 130] . assuming the AdS/CFT 
correspondence, for the strong coupling g — > oo behavior of the scaling function f(g): 

/string^) = 2 V2g - - log(2) + O(-) , (13) 
ii g 

where 2^2g = c.f. ©. The leading 0(g) = 0{yf\) piece is obtained from a 
classical string spinning with a large angular momentum s on the AdS space [2H], while 
the O(g ) = O(A ) term is the first quantum correction obtained in [30J. 

On the other hand, performing the strong coupling limit for our scaling function as 
defined from the integral equation (J8j) with ©,(J9]) (see Section [33]) we do vindicate the 
0(g) asymptotics predicted from string theory: the leading contribution to a is of order 
1/g 2 and eliminates the first term on the r.h.s. of ([7]). However, our analysis is currently 
not precise enough to decide whether or not the sub leading 0(g) term matches (fl3l) . We 
hope to present a more complete solution of the strong coupling problem in future work. 

2 The alert reader will notice that transcendentality could still be preserved if (3 turned out to be a 
rational number times C(3)- The important point is that our calculation leads to a detection mechanism 
for BMN-scaling violation. If a future field theory calculation finds (3 ^ BMN scaling breaks down. 
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2 The Factorized Two-Loop 5 1(2) S-matrix 



In this preliminary chapter we will recall the Bethe ansatz for the sl(2) sector of M = 4 
twist operators (jlj) at one loop p2] and beyond [HI [20] . We will then derive it at two 
loops by Feynman diagram computations, successfully checking part of the conjecture 
of [IB]. A complimentary two-loop approach was recently accomplished by Zwiebel [19j, 
who worked out the full (even non-planar) dilatation operator of the bigger su(l,l|2) 
sector by algebraic means, and also demonstrated the emergence of the two-loop two- 
body S-matrix of [T5] . 

The Bethe ansatz is obtained through the diagonalization of an integrable spin chain, 
whose Hamiltonian is equivalent to the dilatation operator. For a general introduction 
into this technology see [31]. The states of the spin chain are represented by removing 
the trace from the gauge theory states. With s\ + s 2 + ■ ■ ■ + sl-i + sl = s, one has 

"K((V^Z)(V s >Z) . . . (V s ^Z)(V^Z)^j — > \ Sl , s 2 , . . . , Si _ a , s L ) , (14) 

corresponding to a chain of length L. The Sj are the spins of the chain, and can, if s is 
sufficiently large, take on any value due to the noncompact character of the sl(2) sector. 
The Hamiltonian of the chain acts on this state space. Anomalous dimensions A are 
then related to the energies E(g) (i.e. the eigenvalues of the Hamiltonian) through 

A = L + s + g 2 E(g). (15) 

Recall the one-loop Bethe ansatz for sl(2), corresponding to a XXX_i nearest- 
neighbor spin chain where the subscript indicates a non-compact spin — | representation 
of*[(2): 

^4) L = n^^i, *=i,...,.. tie) 

&* 

The cyclicity constraint and the one- loop energy Eq := E(Q) are 

k=l Uk 2 k= l U k + 4 

For a pedagogical derivation of these expressions from the Hamiltonian, using coordinate- 
space Bethe ansatz, see [IB]. A rigorous proof, for any representation of st(2), may be 
found in |33j. 

The conjectured asymptotic all-loop Bethe ansatz for sl(2) [20] is then obtained by 
"deforming" the spectral parameter u, where the deformation parameter is the Yang- 
Mills coupling constant g: 

u± ~ = x ± + 9 ncA 



2 2x ± 



It reads: 



x k 



x 

3 



TT^ X JL 1 ^l 2x \ x i ; k = l,...,s, (19) 



k) | = i X t ~ x j 1 ~9 2 /2x k x] 
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with the new cyclicity constraint and the asymptotic all-loop energy E(g) being given 
by 



x + 



n£ = l and E{g) = Y. \i* ~ i= )■ < 2 <» 



k=l X k k=1 \ x fc x fc 



It generalizes a three-loop Bethe ansatz first proposed in [18J. Very recently, much 
additional support of the ansatz was obtained in [21J. It should be noted, however, that 
we still cannot currently prove that the ansatz (|19|) . (|20p really diagonalizes the gauge 
theory at four loops and beyond. The reason is that we do not know how to fix possible 
"dressing factors" (see [21] and references therein.) 

Directly proving the higher-loop ansatz from the gauge field theory is hard. For all 
loops it will surely require ideas that go far beyond "summing up Feynman diagrams". 
To illustrate the complexity we will nevertheless derive the Bethe ansatz at two loops by 
traditional methods. Actually, we will be able to find the S-matrix, and we will succeed 
in checking two-loop factorization of the three-body problem. This is, according to [32J, 
a strong test for integrability. Completing the proof would require to demonstrate the 
factorization of the s-body problem for arbitrary s, which we have not attempted to do. 



2.1 One-Loop Bethe Ansatz and Three-Body Factorization 

The sl{2) sector contains composite operators built from only one complex scalar field Z 
of the A/" = 4 SYM set of fields and the Yang-Mills covariant derivative D^. The operators 
are taken to carry symmetric traceless irreps of the Lorentz group. We may project all 
indices onto the complex direction z = (x% + ix 2 ) / V2, which guarantees symmetrization 
while the trace terms automatically vanish. 

Single trace operators of this type have a natural description as spin chains: each field 
Z is interpreted as an empty site which may be occupied by any number of derivatives 
D z . The spin chain Hamiltonian 

= J2 W S 0) ( 21 ) 

1=1 

involves a nearest neighbour interaction Tif 1 that cyclically acts on all sites of the chain 
of length L. Alternatively, one may consider the asymptotic case, i.e. an open chain of 
infinite length. The Hamiltonian can transfer derivatives and it is conveniently expressed 
by matrices containing amplitudes for such processes. 

The one-loop Hamiltonian was worked out in [33]: let us denote the number of 
derivatives on two adjacent sites as {si, s 2 }. Then 

nf\{ Sl , s 2 } - { Sl , s 2 }) = h( Sl ) + h(s 2 ), (22) 
^({si, s 2 } -> {si - d, s 2 + d}) = -7^7 

where h(k) is the A;-th harmonic number. The matrix elements refer to a basis in which 
{si, s 2 } is divided by s x ! s 2 \ in order to account for the indistinguishability of the deriva- 
tives at each site. 
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The Bethe ansatz rests on the observation that the derivatives D z behave like particles 
(or "magnons") whose motion is governed by a discrete Schrodinger equation. Let us 
assign a position Xj and a momentum p t to each magnon. One constructs a wave function 
for each magnon number s: 

s = l: J2 X1 * (0 \xi)\xi), (23) 
s = 2: Ex!<x a y { °\x 1 ,x 2 )\x 1 ,x 2 ), 

S = 3 : Ea^a^s ^ (0) (>1> ^2, X3) kl ? ^2, S3) 

Here \xi) denotes a state with a magnon at position x\ etc.. The Hamiltonian reshuffles 
these "kets" as it can shift magnons. On the other hand, the kets form a complete set of 
states whose mutual independence one may use to transform the Schrodinger equation 

H (o) J- * (0) (x u x 2 , ...)\x 1 , x 2 , ...) = E<® * {0) (xi,X2,...)\x h x 2 ,...) 

X1<X2<... Xl<X2<... 

(24) 

into a difference equation on the wave function \l/( ). We find e.g. for only one magnon 
2 ¥°\ Xl ) - ¥°\ Xl - 1) - ¥°\ Xl + 1) = £ (0) ¥ \ Xl ) . (25) 
This can be solved by Fourier transform: 

= e ipixi , £ (0) = 4 sin 2 (^) . (26) 

The one-magnon problem thus defines the dispersion law, i.e. the dependence E(p) of 
the energy on the momentum of the particle. It is an essential assumption of the Bethe 
ansatz that the dispersion law for several magnons is simply a sum over the contributions 
of the individual pseudo-particles given by (12B]) . 

For two magnons the arguments of the wave function x 2 ) should obey x\ < 

x 2 in order to avoid over-counting. Since the one-loop Hamiltonian is a two-site inter- 
action, the plane wave solution remains valid when the separation of the magnons is 
greater or equal two. The corresponding difference equation looks in fact like two copies 
of (J2SD: 

2 ¥°\x u x 2 ) - ¥°\ Xl - 1, x 2 ) - ¥°\ Xl + 1, x 2 ) (27) 
+ 2 ¥ 0) ( Xl , x 2 ) - ^ {0) (x!, x 2 - 1) - ^ (0) (xi, x 2 + 1) 
= g 2 M(E^( Pl ) + E^fa))*^, x 2 ). 

It is a special feature of the Hamiltonian ( 1221) that this equation remains valid when 
X\ = x 2 — 1. However, we do find a new equation when x\ = x 2 [T8] : 

^ ¥°\x u x 2 ) - ¥°\ Xl - 1, x 2 ) - l - *(°)(xx -l,x 2 - 1) (28) 

+ 5 ^(X!, X 2 ) - *<°>(zi, X 2 + 1) - i + 1, X 2 + 1) 

= g 2 M(E ( -°\p 1 ) + ^(pa))*^^!, x 2 ) . 
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A simple plane wave does not obey this equation, but we can solve by an ansatz of the 
form 

¥°\x u x 2 ) = e ipiXl + ip2X2 + S {0) (p 2 , Pl ) e ip2X1+ipiX2 . (29) 

The physical intuition behind the last formula is that the particles may scatter by ex- 
changing their momenta; the second plane wave is related to this, whereby the factor 
is called the scattering matrix. It can be determined from ( 1281) : 

„ipi+ip2 _ Opipi i 1 

S^(p 2 , Pl ) = -f— Z€ . +i (30) 

Note that the two plane waves in ( 129]) (with straight and flipped momenta, respectively) 
are independent as functions. Equation ( 1281) therefore yields two conditions, although 
they are equivalent in this case. 

For three magnons one writes an ansatz involving a wave function *ff(°\xi, x 2 , x 3 ) 
subject to x\ < x 2 < x 3 and proceeds to set up difference equations. As before, the 
magnons do not feel each other when xi + 1 < x 2 < £3 — 1. One might expect special 
behaviour when x\ = x 2 — 1 or x 3 = x 2 + 1, but due to the structure of the Hamiltonian 
this actually does not yield any new conditions. Thus it remains to investigate the cases 

(i) x\ = x 2 < X3 , (ii) xi < x 2 = X3 , (Hi) x\ = x 2 = £3 . (31) 

We write an ansatz which straightforwardly generalizes the two-magnon formula (1291) : 



^ (0) (Xi, X 2 , X 3 ) = e ipix l +ip2X2+ip 3 x 3 +S m e ipix 1+ ip 3 X2 + ip2X 3 ( 32 ) 

+ e i P2X\+ipiX2+ip 3 x 3 _|_ g(0) e ip2Xi+ip s X2 + ipix 3 

21_3 231- 

_|_ e ip 3 xi +ip ± x 2 + ip2X 3 _|_ e ip 3 xi +ip2X2 + ipix 3 

31-2 321- 

If X\ = x 2 < £3, the difference equation can be separated into three independent pieces 
according to which momentum multiplies X3 in the exponentials. The case X3 = x 2 > x\ 
obviously allows for a similar distinction w.r.t. x\. Five of the resulting six equations 
are independent so that one may solve: 

S& = SW(p 3 ,p 2 ), (33) 

SgJ = S^(p 2 , Pl ), 

4 3 } i = S^(p 2lPl )S^(p 3 ,p 1 ), 

Si% = S ( -°\p 3 ,p 1 )SW(p 3 , P2 ), 

= S^(p 2 ,p 1 )S^(p 3 , Pl )S^(p 3 ,p 2 ). 

This solution persists when all three magnons coincide, which is again a non-trivial 
consequence of the structure of the Hamiltonian. We see that the scattering remains 
non-diffractive, i.e. the momenta are unaltered while they may be exchanged between the 
magnons. What is more, the three-particle S matrices factor into two-particle processes. 
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2.2 Bethe Ansatz and Three-Body Factorization at Two Loops 



The original Bethe ansatz described in the last section may be generalized to higher 
orders in perturbation theory [18]. To this end one writes a perturbation expansion of 
all relevant quantities, namely the Hamiltonian, the ingoing wave and the S matrix. The 
central topic of this section is to derive the two-loop correction to the S matrix in the 
sl(2) sector directly from the M = 4 field theory, and to check three-body factorization 
to two loops. 

In the appendices \K\ and [B] we derive the two-loop Hamiltonian for one, two, and three 
magnons from a graph calculation using Af = 2 superfields [35] and the SSDR scheme 
(supersymmetric dimensional reduction) [36] . The supergraph formalism is preferable 
because it minimizes the number of Feynman integrals; for the present purpose the 
M = 2 formulation is superior to M = 1 supergraphs. We end up with a manageable 
calculation involving about twenty graphs. SSDR is the best suited regulator since it 
allows one to treat superfields in a version of dimensional regularization^ 

We attack the problem of calculating quantities with open indices by tensor decompo- 
sition and employ the QCD package Mincer j37j to evaluate the resulting scalar integrals. 
The package uses 4 — 2e dimensional vectors so that we explicitly have to symmetrize 
and take out trace terms. This makes the computer algebra very awkward so that we 
have limited the scope of the present work to low magnon numbers. The method was 
detailed in [38] by one of the authors. We will heavily draw upon this reference in the 
appendicesjj Appendix |A] reviews the renormalization of two-loop two-point functions 
in dimensional regularization. In Appendix [B] we introduce operators Di, D 2 which gen- 
erate the singular part of the one- and two-loop two-point functions and we show that 
the second anomalous dimensions are matrix elements of the combination D 2 — 1/2 D\, 
thus reproducing the two-loop effective vertex given in [13], where the renormalization 
of the dilatation operator in dimensional regularization was first discussed. Finally, the 
Di are constructed from the supergraphs and the two-loop Hamiltonian is worked out 
for one, two, and three magnons. 

In this section we display the Hamiltonian as it arises from D 2 — 1/2 D\ alone. 
One can introduce into it a number of gauge parameters which do not appear in the 
difference equations defining the wave function and the S matrix. This freedom is (more 
than) sufficient to make the Hamiltonian hermitian and to make the sum of all elements 
in each row or column disappear, as was the case for the one-loop dilatation operator. 
In Appendix |B] we also give another set of transfer rules which includes the contribution 
of a term —1/4 [hf-, Hq-], which is needed when the dilatation operator is required to 
reproduce the 0(g 2 ) re-mixing of the sl(2) sector operators. This term cannot be made 
hermitian by the aforementioned gauge transformations and thus from the point of view 
of the Bethe ansatz it is maybe best omitted. It is interesting to note however, that the 
commutator term does not change the S matrix, while it seems to make redundant any 
wave function renormalization in the Bethe picture. 

The disconnected pieces of the two-loop combinatorics do not influence the Hamil- 

3 At the current loop order the scheme cannot lead to ambiguities arising from the e tensor. 
4 The calculation of the two-loop anomalous dimension of the twist two operators [H] uses similar 
techniques; it is more fully automatized but renounces on the use of supergraphs. 
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tonian. The connected two-loop graphs can stretch over three adjacent sites. The basis 
elements below denote the number of covariant derivatives at these three sites; we have 
explicitly indicated a factor 1/ (si\ s 2 \ S3!) with which they were rescaled. 

Spin 1 

basis: {{1,0,0}, {0,1,0}, {0,0,1}} 



( - 5 
* 4 



\ 



I 



Spin 2 

basis: {§{2,0,0}, {1,1,0}, {1,0,1}, |{0, 2, 0}, {0,1,1}, |{0, 0, 2}} 
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Let us now focus on the Bethe ansatz. The spin-chain Hamiltonian up to two loops is 

9ymN 



w = £«i 0) + ^wi 2) +---> 



9 2 = 



8rr 2 



(34) 



and it has energy eigenvalues E = E^ + g 2 E@) + The wave functions of the form 

*( Xl , x 2 , ...) = ¥°\ Xl , x 2 , ■■■) +g 2 tt (2) (zi, x 2 , ...) + •• ■ (35) 
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are contracted on the kets \xi, x 2 , ■ ■ ■) defined in Section ( 12. ip . 

For one magnon we may scale away ty( 2 \ The Schrodinger equation 

^ (2) (i) E * (0) (zi) ki> = £ (2) E * (0) ( xi ) ( 36 ) 

i x\ xi 

leads to the difference condition 

_l#(o)( a . 1 _2) + 2*^(ari-l) -3^ (0) (xi) + 2*^(^ + 1) - 1 *(°) - 2) 
2 2 

= ^(2)^r(0)( Xl ) (37) 

which can again be solved by Fourier transform: 

= e tpix \ £ (2) = -8sin 4 (^). (38) 

2 

Hence the solution of the two-loop one-magnon problem yields the correction to the one- 
loop dispersion law (126|) for the magnon energy E(p) = E^°'(p) + g 2 E^ 2 \p) + 0(g i ). It 
is identical to the one of the su(2) [JS] and su(l|l) [IE] sectors, and consistent with the 
proposed all-loop dispersion law of [TU] : 

E(p) = ±(^l + 8g* sin 2 (|)-l). 

The lowest order of the two magnon problem was discussed in the last section. The 
two-loop part of the Schrodinger equation reads: 

J^nP * m {x u x 2 )\x u x 2 ) (39) 

i x±<X2 

+ J2 H ? ] E * (0) (^,^) \xt, X 2 ) 
i x±<X2 

= (£ (0) (pi) + E^\p 2 )) * (2) (^2)|x 1; x 2 > 

Xl<X2 

+ (E^( Pl ) + E^{p 2 )) * { °\xi, X2) \xi, x 2 ) 

X\<X2 

The resulting difference equations are perhaps not particularly elucidating. We will 

rather comment on how to solve the system: the interaction length of the two-loop 

(2) 

Hamiltonian 7i\ is three. The two magnons must therefore behave as free particles 

(2) 

when x\ < x 2 — 2. Thanks to the special form of Tt\ the same difference equation still 
holds when x\ = x 2 — 2. The cases of interest are thus 

(i) xi = x 2 — 1 , (ii) x\ = x 2 , (40) 

which both lead to new equations. In order to satisfy both conditions we must allow for 
a correction not only to the S matrix but also to the ingoing wave function. Let 

^(Pi,p 2 ) = (1 + 9 2 5 xl>x J(p^p 2 ))^ +ip ™. (41) 
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The wave function renormalization ("fudge factor") is local. We write the ansatz 

*{ Xl , x 2 ) = V(pi, p 2 ) + S(p 2 , pi) ^(pa, Pi) , ^ = ^ (0) + 9 2 S {2) (42) 

whose expansion in the coupling constant defines ty^ 2 '. 

Case (i) gives a condition relating /(pi, p 2 ) to f{p 2 , Pi)- Substituting this into (ii) 
makes the fudge factors disappear from the equation so that we can solve for the S 
matrix: 

8* sin(f ) sin(^) sin(f ) (sin 2 (f ) + sin 2 (f )) 
lP2 ' Pl> ~ (sin(^) + 2i sin(f ) sin(f )) 2 1 1 

This result nicely confirms the conjecture for the two-loop S-matrix of the sl(2) sector 
in [IB]. 

The wave function renormalization / is not fully determined. It is tempting to assume 
it to be symmetric under the exchange of pi, p 2 since the magnons are indistinguishable. 
In this case we find 

ft \ ■ *( Pl \ i • 2/P2N 1 . 2/Pl+P2\ ( A A\ 

f(Pi, P2) = sin {—) + sin {—) - - sin ( — ) . (44) 

The alternative choice for the two-loop Hamiltonian from Appendix [B] yields f(pi, p 2 ) = 
if / is symmetric. 

The discussion of the two-loop three magnon scattering combines elements of the one- 
loop three magnon case with the two magnon situation described in the last paragraph. 
We write for the ingoing wave 

1>(pi,P2, Pa) = (l + g 2 M (5 Xl)X2 l(pi, p 2 , p 3 ) + S X2>X3 r(px, p 2 , p 3 ) (45) 



V>(jPl, P2, P3)) ) 



tpixi + ip2X2 + IP3XZ 



and make the ansatz 



*(ari, x 2 , x 3 ) = V(Pi, P2, Pa) + Si 32 ip(pi, Pa, P2) (46) 

+ S 213 ip(p 2 , pi, p 3 ) + 5*231 ip(p 2 , Pa, Pi) 

+ S 312 1p(p 3 , Pi, p 2 ) + 5321 V(P3, P2, Pi) , 
C _ c(0) , „2 o(2) 

&ijk — ^ijk y ^ijk ■ 

As might be expected by now, the free situation must arise when x\ + 2 < x 2 < x 3 — 2 
but in fact nothing changes when X\ + 2 = x 2 01 x 2 = x 3 — 2. We thus have to discuss 
the cases 

(i) xi + 1 = x 2 < x 3 - 1 , (47) 

(ii) xi = x 2 < x 3 — 1 , 
(Hi) x\ + 1 < x 2 = x 3 — 1 , 
(zu) X\ + 1 < x 2 = x 3 , 

(v) Xi + 1 = x 2 = x 3 - 1 , 
(w) Xi = x 2 — x 3 - 1 , 

(vii) xi + 1 = x 2 = x 3 , 

(viii) Xi = x 2 = x 3 . 
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In the first four cases only one x has disappeared whereby one may use the functional 
independence of the various exponential factors to organize each difference equation into 
three separate constraints. Cases (i) and (ii) are equivalent to a two-magnon problem 
with positions X\ and x 2 : one may solve (i) for three conditions relating l(pi, p 2 , Pz) to 
l(p2, Pi-, pz) etc. and then substitute the three equations into (ii). This eliminates the 
left fudge factor I from the equations. Likewise, we can use (iii) to eliminate the right 

(2) 

fudge factor r from (iv). We are left with six equations on the five S\- k matrices. A 
unique solution exists: 

(Pa, P2) , (48) 

(P2, Pi) , 

(P2, Pi) ^ (2) (P3, Pi) + SW(p 2 , Pi) SMfa, Pi) , 

(p 3 , Pi) S^(p 3 , pa) + S™(p 3 , p x ) 5(°)(p 3 , pa) , 

(pa, pi) 5(°)(p 3 , pi) S^(p 3 , p 2 ) + 5-(°)(p 2 , pi) S (2) (p 3 , px) ^(pa, p 2 ) + 

(p 2 ,Pl)5(°)(p3,p 1 )^ )(p3,p 2 ). 

In other words, the complete S matrix S = + g 2 factors into two-particle 

processes also at two loops. 

Once knowing that 7i- (3) reproduces a two magnon problem when only two ar- 
guments coincide, it is natural to put Z(pi, pa, P3) = /(pi, P2) = r(p 3 , Pi, P2) and so 
on. With these identifications the cases (i) and (iii) reduce to the condition on the two- 
magnon fudge factor / found earlier. Of our remaining cases (v) is empty while the last 
three all lead to one and the same condition on the ultra-local fudge factor u. There is 
not enough information at this loop order to solve for u — again, one may speculate that 
it should be chosen so as to make the ingoing wave symmetric when all three positions 
coincide. The solution is then similar to (PHI) if the two-loop Hamiltonian is as defined 
in this section, or it vanishes for the alternative choice of Tif^ 1 from Appendix [B] 

In conclusion, our analysis confirms the possibility of extending the si (2) sector Bethe 
ansatz to the two-loop level. It proves the functional form of the two-loop S matrix 
conjectured in [IB] , and it shows that the three- magnon S matrix factors into two-particle 
blocks. 

3 The Asymptotic All-Loop Large Spin Limit 
3.1 One-Loop Large Spin Limit 

Consider the one- loop Bethe equations (116l) . (117j) in the large spin limit s — > 00. This 
problem was solved in great detail in the context of Reggeized gluon scattering for the 
very similar case of a noncompact sl(2) spin= representation, i.e. for a XXXo Heisen- 
berg magnet, in [22]. The changes required to treat our present case of noncompact sl(2) 
spin= — I are minor. Here we will proceed in a slightly different fashion as compared to 
|39j . where methods involving the Baxter-Q function are employed. The reason is that 
the higher loop generalization of the Baxter function is not yet known. We will therefore 



S. 



°132 

(2) 
213 

o(2) 
°231 

c(2) 
°312 

o(2) 
°321 



S^ 2 

C(2 
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directly work with the one- loop Bethe equations ( |T6l) .f |T71) . which nicely turn into a (sin- 
gular) integral equation in the large spin limit. Our method will then be extended to the 
asymptotic all-loop equations (jT9"j) . (1201 in the next section. Interestingly, the effective 
higher-loop integral equation will turn out to be non-singular. 

Much intuition may be gained from the fact that the twist L = 2 case is, at one loop, 
explicitly solvable for arbitrary spin s, cf. Appendix O Studying this solution one finds 
that the Bethe roots are all reaO and symmetrically distributed around zero. The root 
distribution density has a peak at the origin (in particular, there is no gap around zero) 
and the outermost roots grow linearly with the spin as max{|iifc|} — > s/2. We therefore 
introduce rescaled variables u, and a density po(u) normalized to one: 

— — — > u with p{u) = - 5 (u — — ) and thus / Po{u) = 1 . (49) 

s tl s J -~ b 

We now take the usual logarithm of the Bethe equations ( JTBT) and multiply either side 
by -i: 

- i L log (^±i) = 2 rrn k - i V log u ^ u i-\ . (50) 

\U k -\) U k -Uj + l 

The integers n k reflect the ambiguity in the branch of the logarithm, and may be inter- 
preted as (bosonic) quantum mode numbers. In the case of twist L = 2 there is only one 
state. Its root distribution is real and symmetric under u <-> — u. All positive (negative) 
roots have mode number n = 1 (n = — 1). In the case of higher twist L > 2 there is 
more than one statfi However, for the lowest state the root distribution is again real 
symmetric with n =sgn(w). Since s is assumed large, and u k = O(s) for nearly all roots, 
we furthermore expand (T5*U|) in 1/u: 

L s 1 

— = 2vrn fc -2^ . (51) 

3=1 



In this large s limit the rescaled Bethe roots condense onto a smooth cut on the interval 
[—b, b] on the real w-axis. We may therefore take a continuum limit of floTl) which yields, 
using (T491) . 

= 27re(u)-2 / du' , (52) 

J-l u-W 

where e{u) =sgn(w). In particular, the dependence on L in ( f51i) drops out: The lowest 
state leads to the same large s root distribution, and therefore energy, for arbitrary finite 
twist L. 



5 It may be shown that, in contrast to the su(2) spin=i Heisenberg magnet, the roots of the sl(2) 
spin=-i Bethe equations are, for all L and s, always real. We thank V. Kazakov and K. Zarembo for a 
discussion of this point. 

6 The reader might find it instructive to consult Table 2 of [TB] , where a complete list of the three- loop 
spectrum of the first few states of the sl(2) sector may be found. 
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The singular integral equation ( 1521) is easily solved by inverting the finite Hilbert 
transform with standard methods. The solution for the rescaled one-loop root density is 
then found to be 



_ . _ . 1 1 + V1-4M 2 2 /, — \ , , 

Po(u) = — log . = — arctanh v 1 — 4tr , (53) 

where we have set the interval boundary to 6 = |, as obtained from the density normal- 
ization condition. The result (153*1) of our procedure agrees with the Baxter-Q approach 
of [33]. 

Our derivation is closely modeled after the discussion of [40J; in particular, we refer 
to appendix C of that article. There the "spinning strings" solutions of ( fl6l) . (flTl) . where 
both s and L are large and of the same order of magnitude O(L) = O(s), were studied. 
The difference is that in this case the l.h.s of ( 15 ip is not negligible. The ensuing potential 
L/u on the l.h.s. of (j5Tl) opens up a gap [—a, a] of the root distribution in the vicinity 
of u — 0. The resulting density for the lowest state therefore has compact support on 
two cuts [—b, —a] and [a, b] and is expressible through an elliptic integral of the third 
kind (see eq. (C.8) in [40J). One easily checks that when L — > the gap disappears, 
i.e. a — > 0, and the elliptic density, after rescaling the roots in [ID] by u — > fit in order 
to adapt conventions, simplifies to the expression (1531) . with b — > ~. 

However, the one-loop anomalous dimension as obtained in [HJ] does not reproduce 
the expected logarithmic scaling of (j5J upon taking the limit s/L — > oo. Instead, it be- 
haves like ~ log 2 (s), cf. (E.l) of [10]. This is a classic order-of-limits problem. Assuming 
s, L large with s/L finite, and subsequently taking s/L — > oo does not yield the same 
result as taking s large while keeping L either finite or, at least, L s. For a very 
recent, quite extensive discussion of this fact see [27J. For a recent study of some of the 
fine-structure of the spinning strings limit see [4~T] . 

The correct result is obtained by a careful derivation of the expression for the energy 
in the continuum limit s — > oo. From the right equation in (1171) we find, using 

1 [* j- Po( u 



Eo = - du \ . (54) 



S / i u 2 + 

2 



1.S 



Therefore, as opposed to the limit of |40j (see the expression in (C.4)) it is nonsensical 
to use the unregulated expectation value f dupo{u)/u 2 for the energy. The correct 
expression (1341) is actually related to the resolvent G(u), which is defined for arbitrary 
complex values of u barring the interval [—1/2, 1/2] (this integral is not of principal part 
type) as 

G(u)= T 2 du'^1, (55) 
.1-1 w — u 



through 



Note that this further distinguishes the large spin limit from the "spinning strings" limit, 
where the resolvent generates the full set commuting charges jl2]. One then finds from 
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531 that 



G(u) = z log —=== — - . 57 
VI - 4m 2 - 1 

Using now ( 1561) and taking s->oowe find 

£ = 4 log(s) + O(s ) , (58) 

which is the well-known correct result, as may also be checked directly from the exact 
finite s result E = 4h(s), see ( 11521) . 

3.2 Asymptotic All-Loop Large Spin Limit 

Let us now generalize the analysis of the previous section to the higher loop case. We 
would therefore like to compute the corrections to the one-loop density (1531 and energy 
(15H1) as generated by the deformed Bethe equations (fTUl) . (12U|) . Compelling arguments for 
its validity to three loops were presented in [18] (in particular the equations reproduce 
the conjecture of [ID] based on the QCD calculation [9], and they agree with [2]). Their 
all- loop form was conjectured in [2D]. See also [2T] . 

We begin by rewriting the asymptotic all-loop Bethe equations (1T9|) with the help of 
( 1181) in the following fashion: 

^ + j \ L ( l+9 2 /2(x- k r L 




Let us again take a logarithm on both sides of the equations, and multiply by i: 

2L arctan(2 u k ) + i L log I + = 2nh k -2 } arctan (u k - uj] 

\l + g 2 /2(x+) 2 J .^ /2 



+ 2i 



i E l0 S 2 + ( 6 °) 



j=-s/2 

Here we have also relabeled the s roots u k such that the index k runs over the set 
k = ±1, ±2, . . . , ±|. We have furthermore chosen, for convenience, to employ a different 
choice for the branches of the logarithms as compared to ( 1501) . Whereas in ( 150|) the 
branchcuts run through u k = and u k = Uj, in our alternative choice in ( l60l) the arctan 
functions are analytic at u k = and = Wj. This replaces the "bosonic" mode numbers 
n k of (1501) by "fermionic" mode numbers h k . For the lowest state (the only one for L = 2) 
we have, for even s, 

h k = k + for fc = ±1, ±2, . . . , ±| . (61) 
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To avoid confusion: We are still focusing on the same states, and just chose to change 
the description. 

Let us now proceed in close similarity to the computation of the thermodynamic 
antiferromagnetic ground state of the Heisenberg magnet (see e.g. [S3])- In order to have 
a uniform spacing between the indices of all roots it is convenient to define k' = k — e(k)/2 
such that 

L 2 13 1 

n k = k' + —^-e(k) for k' = ±-, ±~, . . . , ±-^- . (62) 

As s —y oo we introduce a smooth continuum variable x — — . The excitation density 
may now be defined as p(u) = We divide fIBll by s, use (J6"2"]l . replace the sums by 
integrals, and, finally, take a derivative w.r.t. u. Note that we do not rescale u by 1/s. 
Then ( 160]) becomes 

L 1 iL d , fl + g 2 /2(x-(u)) 2 

log ' 



su 2 + \ s du *\l + g 2 /2{x+(u)) 



^ 2.pH + -(L-2) W -2| fc ^ (M _ M/)2 + i 

o- P j i / rf i A-^ 2 /2x+(u)x-K)\ 
+ 2z| 6 ^V(,0^1og( 1 _^_^ + ^ j • (63) 

It is convenient to split the density p{u) into a one-loop piece po(u) and a higher-loop 
piece a{u): p{u) = po(u) + g 2 a{u). Let us, accordingly, also split off from (l6Uj) the 
one-loop contribution 

L -J-, = 2 ., oM + (L - 2) W - 2 /' , (64) 

while the higher (two and beyond) loop part of (I63p becomes 

= 27ra(u) 

1 , , a(u') 
du 



>3C 



(u-u'V + l 



2i n , , d (I - g 2 /2x + {u)x-{u') 



[ du poiu') — log 



g 2 /2x~ 



+ 2,I J M log f'-^'f:^ ) • (65) 
' ; du & \1 - g 2 /2x-(u)x+(u')/ 



We have dropped the second term on the l.h.s. of ( 1631) . as it is easily seen to be suppressed 
to leading order in the large s limit. This reflects the independence of the large s scaling 
behavior of the lowest state on the twist L even beyond the one-loop approximation, as 
long as L C s. We have also extended the range of integration of the second and fourth 
integral in (J6"5"]l from ±s/2 to ±oo, to be justified below. 

As a consistency check of our procedure let us rederive the one-loop solution of the 
previous section from ( |64j) . There we used rescaled variables u = -, and a rescaled 
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density po(u) = s po(u) such that dupo(u) = dupo(u). Using the large s expansions 



2s^ + _I_ 



s (u _ u'Y + jf 



ir6(u)+0(-) 



71 5 (u — u) + 



s \u — u 



V 1 
-7W + ( ^)> 



(66) 
(67) 



where V indicates a principal part, we find from ([6 



= Ati5(u) + 2 j du' 



PoK) 
u — u') 2 



(68) 



which is, since e'(u 



2 5(u), precisely the derivative of the one- loop singular integral 
equation ( 1521) . Note that the L dependence has indeed again dropped out. We therefore 
find the same one- loop result as in the previous section. It should be stressed that, even 
though the kernel in is of difference form, and the interval boundary values tend to 
±oo, it is incorrect to solve this equation by naive Fourier techniques. 

Luckily, however, applying a Fourier transform leads to progress with the higher- 
loop equation (1651) . The reason is that the higher loop density fluctuations cr(u) are 
concentrated in the vicinity of u — 0, i.e. &(u) ^ iff \u\ s/2. This may be verified 
for twist L = 2 operators by using the exact one-loop solution of appendix and 
numerically solving the linear problem of computing the higher-loop corrections to the 
roots of the Hahn polynomials from the Bethe equations ffTUl) . We were thus indeed 
entitled to replace the integral boundaries ±s/2 by ±oo in the second and fourth term 
on the r.h.s. of ( 1651) . The "scale" of the fluctuations <t(m) is set by the third term on the 
r.h.s. of ( 1651) . Let us calculate it, using po(u) = po(u)/s, with po(u) given by ( |53i) : 



2i 

2 

n . oo 
y r=l 

1 d 



d 

du p (u) — log 



g 2 /2x- 



[U)X I u 



s 



2du 



X [ U j 



+ 



du' po(u' 
1 



g 2 /2x~ 
d 



X~[U) 



du 
+ ... , 



U)X^ 

1 



u 



x + (u) r x~(u') r x~ (u) r x + (u') 1 



(69) 



where we have only kept the leading contribution. Note that only the first, r = 1 term 
in the expansion of the logarithm contributes to this result, and we have used, cf. fl55l) . 
(1561) . the relation 



du' po(u') 



x ± (u') 



2 1 

du'p (u) -—j + 

2 ^ 



s \2s 2s ° 



(70) 

which is valid to leading order at large s. It is now clear from (lo^j) that E /s ~ 
4 log(s)/s sets the scale of the density fluctuation a{u) in (IBTj) . We therefore define 
a(u) = -{Eq/s) a(u), i.e. 

p(u) = p (u) - g 2 — a(u) . (71) 
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To this leading order, the density fluctuation does not change the density normalization 
J b b dup(u) = 1, i.e. J^_ b du p (u) = 1 + ... since lim E /s = 0, see (1551) . Then (1651) 

becomes 







2 Ti a(u) 

2 r m, 



1 d 

2du 



1 1 

TV"^ + T 



X [ U ) X~ M 



+2 z / du' <j(u') — — lot 
ou 



g 2 /2x + {u)x-(u') 



1- g 2 /2x~(u)x+(u') 
We now introduce the Fourier transform a(t) of the fluctuation density <j(u) 



(72) 



a{t) = e"2 



due a(u 



(73) 



where we have also included a factor e 2 for notational convenience. Fourier transforming 
e~2 due ltu x equation ( 1721) we find, after some calculation (see appendix IDl). 



= 2tt &(t) - 2tt e~* <r(t) - 2tt e" 



Ji(V2gt) 
V2g 



4tt ^ t e" 1 / dtf K{V2 gt,V2g t') a(t') , 
'0 

(74) 



where the four terms in (1741) correspond, respectively, to the four terms in (1721) . and the 
kernel K is given in terms of Bessel functions by 



K(V2gt,y/2gt') 



1 JxiV^g^MV^gt') - J Q {y/2gt)J x {y/2gi!) 



V2g 



t - V 



(75) 



Note that the Fourier transform only diagonalizes the "main" scattering term in (172"]) . 
i.e. the kernel l/((w — u') 2 + 1). So we are still left with an integral equation. However, 
the higher- loop equation (1741) is, in view of (1751) . and in contradistinction to the one- 
loop equation ( |52l) . non- singular. It may be rewritten in the form (j5J) stated in the 
introduction. Finally, the all-loop energy is found from (1201) to be 



E(g) 



du p(u) 



E -g 2 E 



x + (u) x~{u) 

dua(u) ( — l — - — ^-r 
\x+{u) x (it) 



(76) 



to leading order in s. In terms of the Fourier transformed density <x(t), cf. ( 1731) . this 
becomes (see again appendix iDl) 



E{g) = E (^-Ag 2 



dt*(t) J «f 9t U ^ 



(77) 
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with Eq = 4 log(s) + . . . . Notice that, in line with general expectations, we have just 
shown that the Bethe ansatz of [20] indeed leads to the logarithmic scaling behavior ((Sj) 
to all orders in perturbation theory, in agreement with general expectations, see e.g. the 
discussions in [TU] , [39J , [27J . In view of (15]). ([15]) , (1581) this indeed yields our proposed 
conjecture for the all-loop scaling function f(g) announced in (J7j). The proposed scaling 
function as found from the Bethe ansatz possesses further remarkable properties, to 
which we will now turn our attention. 

3.3 Weak-Coupling Expansion and Transcendentality 

The Fredholm form of the higher-loop integral equation (j8J or ( 1741) is ideally suited for 
the explicit perturbative expansion of the scaling function f(g) of (J5J to high orders. 
Both the inhomogeneous, first term as well as the kernel of (j8J) have a regular expansion 
in even powers of g around g = 0. We may therefore also expand the transformed density 
a(t) in even powers of g and solve (J8J iteratively 




where we have used the following representation of the Riemann zeta function: 

1 f°° dtt n 

c( " +1) = di„ — • (79 > 

Furthermore, the expression for the scaling function ((7j) may also be expanded in a Taylor 
series in g 2 : 

9 i f°° dtt a( r° dtt 3 t >/ , f°° dtt \ 
/( 9 )=4 9 >-Vy o — +g >(2j a — +H(2 )J a _)+.... (80) 

We again use (1791) and we find to e.g. six-loop order 

f(g) = Ag 2 -4C(2)s 4 + (4C(2) 2 + 12C(4))$ 8 (81) 

- (4 C(2) 3 + 24 C(2)C(4) - 4 C(3) 2 + 50 C(6)) g 8 

+ (4C(2) 4 + 36C(2) 2 C(4) - 8C(2)C(3) 2 + 100C(2)C(6) 

- 40 C(3)C(5) + 39 C(4) 2 + 245 C(8)) g 10 

- (4 C(2) 5 + 48 C(2) 3 C(4) - 12 C(2) 2 C(3) 2 + 150 C(2) 2 ((6) 

- 80 C(2)C(3)C(5) + 114 C(2)C(4) 2 + 490 C(2)C(8) - 18 C(3) 2 C(4) 

- 210 C(3)C(7) + 345 C(4)C(6) - 102 C(5) 2 + 1323 C(10)) g 12 + . . . 

It is easy to go to much higher orders if desired (we have expanded to 20-loop order g m ). 
It is seen that the £-loop 0(g 2e ) contribution to the anomalous dimension is a sum of 
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products of zeta functions. What is more, the arguments of the zeta functions of each 
product always add up to the number 2 £ — 2. This is a test of the "trans cendentality 
principle" of Kotikov, Lipatov, Onishchenko and Velizhanin as spelled out in [7J El fTUj . 
and we see that our Bethe ansatz is consistent with this principlqj- Finally it is also seen 
that the numerical coefficients in front of each zeta function product are integers^. 

Note that the expansion ( 1HTT) may be written more compactly when expressing the 
zeta functions of even arguments through powers of tt times rational numbers: 

f(9) = 4/ (82) 
9 

^2 4 

71 q 

3 y 

45 y 
73_ 

,630 

+ (^- 8 -^ 2 C(3) 2 -40C(3)C(5)), 10 

- (43?883 _ 8_ ^ 2 _ 40 
V3742200 15 w 3 v ; v ; 

-210C(3)C(7)-102C(5) 2 ) £ 12 
+ .... 

Clearly each 7r contributes one "unit" of transcendentality. This however obscures the 
integer nature of the numerical coefficients (c.f. footnote 8). 

It is instructive to investigate whether the (BMN scaling-preserving) "AFS" dressing 
factor pi IH ED] for the (approximate, see [4H H51 116]) string Bethe ansatz (fT9l) . (1201 is 
compatible with the transcendentality principle. Possible (BMN scaling-violating) gauge 
dressing factors are briefly treated in the next section 13.41 

The AFS ansatz leads to a modification, at three loops and beyond, of the integral 
equation ([H]) 

a(t) = - r - r K'(y/2gt, 0) - 2g 2 / dt' K'(y/2gt, yftgt) ,a(t') , (83) 
e 1 - 1 L J J 



7 To be more precise, here we have tested a weaker form of the transcendentality principle of [10] . 
The stronger form applies to the finite s case, and states that the indices of certain harmonic sums add 
up to 2 1 — 1. We suspect that our all- loop Bethe ansatz is also consistent with the stronger version, see 
also [IS]. Our finding certainly supports this, as the weaker principle is a consequence of the stronger 
one. It would be exciting to fully prove the latter from the L = 2 finite s Bethe equations (jTHJ) , (J^TJJ) . 

8 Actually, with our convention (J2|), higher terms beyond the order we have printed in (|81|) develop 
powers of 2 in the denominator. We however checked up to order g i0 that our scheme yields indeed integer 
numbers in front of the zeta-functions if g is rescaled as g — > V2g, which is Lipatov's et.al. convention. 

9 Our motivation here is not so much string theory as such (in particular we investigate the dressing 
factor at weak coupling, while its original design demands strong coupling) but rather the fact that 
this type of dressing factors are known to naturally appear in certain variant, asymptotically integrable 
spin chains [Uj. While these studies were done for compact magnets, it is likely that they may be 
generalized to the non-compact case of interest in this paper. The variant models tend to violate the 
Feynman rules of the gauge field theory, which is our main motivation for investigating whether they 
preserve the transcendentality principle. 
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where the modified kernel K', see appendix O, reads 

K'(y/2gt, V2gt') = K{V2gt, V2gt') + V2gK(V2gt, V2gt'), (84) 

with 



km = «Mt)Mt)-Mt)Mt)) k(ti0) = m. (85 ) 

t t u 

The dressing factor then modifies the scaling function f(g) — > f(g) + 5f(g) in the fol- 
lowing fashion: 



5f(g) = Oxg 2 (86) 
-4C(3)£ 6 

+ (^vr 2 C(3) + 20C(5))/ 

- (| vr 4 C(3) + f vr 2 C(5) + 105 C(7) - 4 C(3) 2 ) g w 

/ 71 7Q 
+ (3T5 K C(3) + 30 * C(5) + ^ C(7) " 8 C(3)3 + 588 C(9) 

-2 7r 2 C(3) 2 -36C(3)C(5))^ 12 
+ . . . 

We see that the integrable modification of the long-range Bethe ansatz of [20] by the 
"stringy" AFS [43] dressing factor violates the transcendentality principled, as now the 
arguments of the Riemann zeta functions no longer add up to 2£ — 2. 



3.4 Breakdown of BMN scaling and the Scaling Function 

Here we will demonstrate interesting connections between BMN scaling [28] on the one 
hand and our Bethe ansatz method for the scaling function on the other. It is by 
now rather firmly established that BMN scaling in perturbative gauge theory can only 
break down, at four loops or beyond, through a dressing factor of the general type just 
discussed, see in particular [IT], [21]. This happens in e.g. the plane- wave matrix model, 
see [IE]. 

Let us sketch the quantitative derivation of this effect, restricting ourselves for sim- 
plicity to four loops, where its detection might still be within reasonable reach of sophis- 
ticated field theory methods, maybe along the lines of [2] . 

The first modification of the asymptotic Bethe equations of [20] which is still consis- 
tent with current knowledge on the integrable structure of Af = 4 gauge theory would 

10 For the gauge theory ansatz the transcendentality principle is a consequence of scaling: the argu- 
ments of potential and kernel in (|74|) are \plg t, \/2g t' so that the order in g is linked to the total power 
of t and t' which defines the level of transcendentality. The string theory ansatz (|83|) breaks the pattern 
only because of the presence of the extra \[2g in front of K' in equation ([84| . Initially this introduces 
a mismatch by one unit; by iteration the effect fans out higher up in the perturbative expansion. 
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lead to the following correction of the higher loop Bethe equations (JT9 



Xl ) L = II % — ~^ \ 9 -^H *\u k , Uj ) (87) 
with 

a 2 (u k ,Uj) = e i/3 » 6 (" 2 ^)93K)-93K)92(« J ))+... ? ( g8 ) 

see (HI [20J, 07], 05], [21] for details^, and the definition of the charges q r (u). The dots 
indicate further terms which might affect five loops and higher. 

A non-zero value for (3 leads to "soft-breaking" of BMN scaling: The two-excitation 
problem can be solved exactly [TT], because the momentum constraint implies 112 = 
—u\. For spin chain length L = J, where J is the BMN R-charge, the different states 
are distinguished by the Bethe roots M l n = ~ cot (jfj)- The higher order corrections 
similarly come out in terms of trigonometric functions. The string spectrum A — J is 
reproduced by Taylor expanding in 1/J when n is small: 

A-J = 2+8g 2 — -lQg 4 — +64g 6 — -320/ — -512/^ — +. . . 



\ J J \ J J \ J J * V J J * J V J 

(89) 

The dots stand for higher orders in g 2 and, at any given order, terms subleading in 1/J. 
We see the emergence of the effective coupling constant g 2 / J 2 in the first four terms of 
the last formula, while the last term has g 8 / J 7 , so that it diverges in the BMN limit 
g, J — ► 00 with g/J fixed. 

The modified Bethe ansatz ( 1H71) requires replacing the kernel in ( 1751) by 



k(V2g \t\, V2g \t'\) - K(V2g \t\, V~2g \t'\) +2(3(V2g) U^2 g\t\)UV2 g\t'\) + 

(90) 



The third term on the r.h.s. of ( ff2l) becomes 



ld_\ 
2du) 



1 1 
rv^ + — 7 



X^[U) X [u 
or, after Fourier transforming (c.f. third term in ([73 



00 4 ^fc(«) +■■■, (91) 



_ 2 ^-' 1 MVigt) +2 f3g2j 2 (V2g t ) \ . .. mi 



V2g 

This modifies the four-loop 0(g 8 ) term of the scaling function ( l82l) to 

^tt 6 -4C(3) 2 + 8/3C(3)) • (93) 



11 The detailed argumentation which allows to draw this conclusion is actually rather subtle and 
requires putting together various results. The main steps are: (1) The three- loop Bethe ansatz is solidly 
known. (2) The structure of the four-loop Bethe ansatz is also known, up to the term involving [3 in 
([87]) . in the su(2) sector [47]. (3) The multiplicative modification affecting the su(2) sector as in ([87]) 
must also multiplicatively affect in the same fashion the sl(2) sector, as first conjectured in [TB] and 
later proved in [21] . 
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Figure 1: The density of Bethe roots at weak (left) and strong coupling (right). 



Note that transcendentality is violated unless ft is a rational number times ((3) (or 7r 3 ). 
A particularly curious case would be ft = |C(3), which would lead to the much simpler 
four- loop answer — ^ 7r 6 . 

Note that such a modified Bethe ansatz would also change the anomalous dimensions 
of all operators in other sectors. E.g. in the su(2) sector we would find for the length 
L = 5 operator Tr X 2 Z 3 + . . . (this case is actually equivalent to the sl(2) twist three 
operator Tr D 2 Z 3 + . . .) to four loops 

E(g) =A-6g 2 + 17 g A - + 8 ft) g" + . . . . (94) 

It would be very interesting if the modification were non-rationaf^l. Incidentally, we see 
that ft would also rule out the Hubbard Hamiltonian as a candidate for the su(2) 
dilatation operator beyond three-loop order, c.f. eq.(68) in |22j . 



3.5 Strong-Coupling Expansion and String Theory 

The Fourier-transformed integral equation (1741) does not lend itself to strong coupling 
analysis due to the oscillatory nature of the kernel flTo"]) . We rather return to the config- 
uration space integral equation (172"]) . 

The two diagrams in Figure 1 give a series of plots of the root density for progressively 
higher values of the coupling constant. The left picture shows the weak coupling regime; 
the graphs depict the root density at \/2g = 0, 1/4, 1/2, 1, respectively. The \/2g = 
distribution is the tallest peak. It is given by the Fourier back-transform of the first term 
in (US]): 

ao(n) = I J, , (95) 
4 cosh {ITU) 

All other curves are numerical solutions of (1721) . Augmenting the coupling constant 
makes the peak around u = become wider and flatter. 



12 Clearly the so far proposed Bethe ansatze [5D] also lead to a transcendentality principle at weak cou- 
pling: If we assign, in accordance with the meaning of the word, transcendentality degree zero to rational 
or algebraic numbers, then weak coupling dimensions of operators carrying finite charges (i.e. without 
taking limits of large R-charges or large spin quantum numbers) are always of zero degree in the cur- 
rently proposed ansatze. On the other hand, zeta functions do appear naturally in individual higher 
loop Feynman diagrams, and, from this point of view, might well appear in high order contributions to 
anomalous dimensions. 
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In the second diagram we plotted 2g 2 a(u/(V2g)) for V2g = 1, 4, 16, 64. With 
increasing coupling the graphs rise; they develop peaks at ±1 while the middle parts 
tend to l/ir. On undoing the scaling we would nevertheless recover the tendency seen 
at weak coupling, i.e. the support of the root density roughly stretches to the interval 
[— y/2g, V2g] within which the density tends to 

"-("» = dp- (96) 

Note that the constant function u(u) = l/(2irg 2 ) is an exact solution of (172"]) if the 
support is extended to the entire real axis (likewise a(t) = 5(t)/g 2 is a solution of (174]) ). 
Furthermore, a^iu) would exactly cancel the leading 0(g 2 ) contribution to the scaling 
function ( 1771) . thus yielding the 0(g) asymptotics expected from string theory. 

Numerically, we could confirm the cancellation of the 0(g 2 ) part of the scaling func- 
tion up to an error of a few per cent, but reliable predictions for subleading terms 
remained out of reach. It is indispensable to understand the strong coupling regime by 
analytic means. We hope to clarify the issue in future work. 
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A Two-point Functions in the si (2) Sector 

A.l Perturbative CFT in the Dimensional Reduction Scheme 

We shall restrict our attention to leading N (planar) two-point functions of single trace 
operators in the $1(2) sector. For any given spin chain with length=twist L there are 
many distinct operators differing in the total number of derivatives and their positioning 
on the sites of the chain. 

Renormalization must be done in such a way as to un-mix these states and to make 
their correlators finite. The theory is then seen to be conformally invariant; for example 
the two-point function of a renormalized primary operator of spin s has the form [49] 



(P S (1)P S (2)) = gCgS^f^^jW^ , (97) 

where 

J M „(x)=V,-2^ (98) 
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is the inversion tensor, and the fi and v indices are separately made traceless and sym- 
metric. Knowledge of any one term in the product of inversion tensors is sufficient to 
reconstruct the full correlator. In [3B] we considered the term with no 77 symbol, because 
we were interested in a minimal set of graphs (trace terms are potentially more divergent 
and there are also a few Feynman diagrams which always carry at least one power of r\ on 
dimensional grounds). In the present work we wish to construct the asymptotic two- loop 
dilatation operator in the si (2) sector. The task is greatly simplified by focusing on the 
pure trace terms, because these obviously cannot exist between operators of different 
spin. This property becomes important when subtracting out disconnected parts. 

As before, we use Af = 2 superfields and regularize by SSDR (supersymmetric di- 
mensional reduction) [36] in x-space. This amounts to doing the superalgebra as in four 
dimensions, while the underlying scalar propagator is modified as in standard dimen- 
sional regularization: 

(Z(1)Z(2)) = ^-{nx\ 2 y, co = ~, ni<Z(l)Z(2)> = 5(z a ), (99) 

although we suppress the mass scale \i throughout the article^ The tree-level correlators 
of operators of length L and spin s thus contain the x-space structure 

X(L,s) = . - ' z * n r- (100) 

v ' (-47r 2 ) L (xf 2 ) i ( 1 - e )+ s v ; 

and a whole series of terms with x\ 2 with open indices, which may be recovered by 
appealing to conformal invariance. 

In order to extract the one- and two-loop anomalous dimensions we must keep track 
of the leading and sub-leading order in the e expansion of the bare correlators: 

( Oi Oj > = X{L, s) [ (T 0ij + eT Uj ) + g 2 (A 1Uj - + A mj ) (x 2 12 ) e 

+ 9 4 (A 2 2ij ^ + A 21ij - + A 20ij ) {x\ 2 ) 2e + . . . ] (101) 

where the Yang-Mills coupling constant is dressed byM 

and the fractional powers of x\ 2 arise from the integration measure in the Feynman 
graphs defining the one- and two-loop contributions. 

Consistency of M = 4 as a conformal field theory grants that T , An, A 22 are simul- 
taneously diagonalizable. In a diagonal basis {Oi} they obey 

An = Ti T , A 22 = l -Y\T Q . (103) 



13 Our discussion of the renormalization of conformal correlators in x-space using the SSDR scheme 
is built upon the works [T7 | [50 ] [38] . 

1 We deviate from the convention in [38j by a coupling constant rescaling so as to be more in line 
with the literature. 
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Here Ti is also diagonal and contains the one-loop anomalous dimensions 71 j. 
The divergences are removed by introducing Z matrices of the form 

Z = R + g 2 B + g 4 (d - + Co) + ... (104) 

e 

where R is diagonal and has as its entries the Z-factors for the individual operators 

Zi = +g "27 + 9 ^ + IT) + ••• (105) 

while 5, G have zero on the diagonal. The Z factors and the anomalous dimensions are 
determined from the bare two-point functions by imposing 

F = Z{O0)Z ] (106) 

where F is again diagonal and is defined by the renormalized two-point functions 

fi = X(L,s)\ t=0 (a 0i + g 2 a u + g A a 2i ) (a^)-* 2 ^-* 4 ^ + ... (107) 

To be more precise, we demand that both sides be equal at each order in g 2 up to positive 
powers of e. The resulting system of equations does not completely fix Gi, G , so that 
we limit our scope to the determination of R, B, a j, aii, ju, j2i- We may thus drop the 
constant part A20 of the g 4 two-point functions from our analysis. 



A. 2 Graphs 

We exploit the M = 2 superfield formalism in order to minimize the number of Feynman 
diagrams. For a quick review of the essentials of the formalism and expressions for the 
graphs we would like to refer the reader to [38J, where two-loop two-point functions 
of operators of length three are discussed. Our notations and conventions are in fact 
borrowed from that work; in particular, the article contains a list of graphs upon which 
we draw here. However, in [38J the (x z x z ) s term of the two point functions was used, so 
that some graphs could be omitted because they always come with r\ zz . 

At order g 2 , we additionally have to take into account a graph F (see Figure 3 below) 
in which a free vector line goes from the connection in D z on the left end of the two-point 
function to that in D z on the right (the Feynman gauge vector propagator is proportional 
to 77). Correspondingly, there is an 0(g 4 ) graph consisting of the same free line paired 
with the divergent one- loop graph Gq. On the other hand, we do not need to consider 
the combination of the free vector line with the "BPS-like" 0(g 2 ) integral Bq since this 
configuration stays finite. Next, in [38] we could drop the product Go * B as Go only 
has a simple pole (in x-space) while the part of B Q without 77 is a contact term also 
when there are partial derivatives on the outer legs, i.e. it is always 0(e). Terms in B 
which involve 7] zz are finite, i.e. 0(1), so that in the present context the product Go * B 
becomes relevant. 

With respect to the genuine two-loop integrals there are not many changes: the 
finiteness of some terms which we dropped from graph G3 remains guaranteed and hence 
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we may take over the simplified sum G% + G4 given in formula (61) in [3E]. The "BPS- 
like" graphs behave in the same manner as B : the part without rj is a contact term and 
the other parts are finite. They can still safely be omitted. 

A first difference is that graphs Gio and G n start to contribute: before, the poles 
from these graphs cancelled in the sum over all diagrams within each class; this is not 
the case in the new situation 1^1 

But there are also three genuinely new graphs: 




Figure 2. Additional graphs at order g 4 . 



Like in the pictures in [38J we have omitted free matter propagators. Point 1 is on 
the left and point 2 on the right of the graphs. The connection carries the indices [l and 
v there, respectively, while the connection at 1' has index p. The lines are split only for 
convenience of drawing — the notation 1' in G\% does not refer to a new point. It was 
introduced in order to distinguish the two vector propagators joining the cubic vertex 
from the left. 

After the evaluation of Grassmann- and SU(2) integrations we find 

G w = (12) [-77^/4], (109) 
G 17 = 77^/2, (110) 
G 18 = 1 [{d x - d v ) u rjj4 - (£>! - ^ 2 ) P W 4 + ( d i' ~ W 4 ] • (HI) 

In the same way as graphs G 3 , G 4 in [38] occur together, we may add Gi 6 and its mirror 
image G\q into G\\ because their combinatorics is equal: 

G n + G 16 + G 1G = (12) [-3,13^3 + ^u(d 13 -d 23 ).(d u -d 24 )/A (112) 

+ 77^(n 14 + n 24 )/4 + ...] 

(The dots indicate omitted finite terms.) 

The rest of the calculation proceeds along the same lines as before (appropriately 
adapted to the new tensor component), i.e. the reconstruction of the Fourier transform 
of integrals with open indices from projections with the total momentum q and the 77 
symbol, which are built in Mathematica and evaluated by the Mincer package [37J. 

15 We point out an error in formula (55) in the original version of [38) : two parts of the integral were 
added with a wrong relative sign. The cancellation of the associated poles in the calculation of [3H] can 
be verified for both "halves" on their own so that the mistake did not show. The correct expression for 
Gio is: 

do = (12) [^14^23 +^n 34 /4 - (1 - (12)(l-2-))0 Ml3 ft,24] (108) 
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A. 3 The Length 3 Spin 3 Mixing Problem 



As an illustration of what has been said before we re-examine the mixing of the length 
three operators 

K s 2 , s 3 } = Tr((D s /Z)(D?Z)(D s /Z)) (113) 
at leading order in N. In particular, the spin three mixing problem involves the operators 



B = {{3,0,0}, {2,1,0}, {1,2,0}, {1,1,1}}. 



(114) 



The one-loop logarithms and the constant order T of the tree-level correlators are diag- 
onalized by choosing the directions 

O = {1,3,3,2}, (115) 
JC = {1,-1,-1,-2}, 
Vi = {2, -9, -9, 24} , 
V 2 = {0,1,-1,0} 

relative to the basis B. Note that Vi, V2 have identical first anomalous dimension and 
therefore the eigenspace may be spanned by any two independent directions. We have 
split into an even and an odd part under reversal of the trace; as a consequence V2 
decouples from the other operators. Renormalization in the MS scheme outlined above 
yields the anomalous dimensions 
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up to terms of 0(g 6 ). The individual Z{ are given by the anomalous dimensions in 
the standard way. As explained above, the system does not entirely determine the C 
matrices, while we can fix B: 
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(117) 



The parameter a is not calculable from our system of equations because the anomalous 
dimensions of V\, V2 are degenerate. We may put it to zero bearing in mind that an 
arbitrary remixing of the two operators is possible. 

The anomalous dimensions and the entries of B are independent of whether we cal- 
culate the (r/zf) 3 terms as outlined in this article or the (x z Xz) 3 part of the correlators 

ii| 



as in [38], although now in MS 



16 In [3S] we deviated from the strict MS prescription by choosing the Basis B in an e dependent way, 
so that Ti became diagonal as well. This had the advantage of decoupling the protected operator O. 
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The operators Vi and V 2 are conformal primaries of spin three. The numerator of 
their renormalized two-point functions should contain three powers of the inversion tensor 
Jzz — Vzz — ^x z Xz/x 2 , and correspondingly we find that the normalization of the (x z x z ) 3 
terms differs by —8 from that of the {f] Z z) 3 part. The operator /C is a first derivative of the 
primary JCq [38J. The normalizations of the two terms in (KK) are indeed consistent with 
being derivatives of a common spin two two-point function; similarly for the protected 
operator O = 1/3 D\ {0, 0, 0}. In conclusion, in this example renormalization works in 
the same way for these two components of the tensor structure. Conformal invariance is 
manifest . 

B The Dilatation Operator and Renormalization 

B.l Matrix Elements of the Dilatation Operator in Dimen- 
sional Regularization 

Suppose there are linear operators L>i, D 2 

DiOi = (^D 1Uj + D 10ij )Oj, (118) 

D 2 3 = (^D 22ij + l -D 21lj )0 3} (119) 

such that 

(O i O j ) g 2 = ({DtO^O^go = (-D u T + (D W T + D X1 T X )) (120) 

(OiO^ = ((D 2 Oi) Oj) g o = (^D 22 T + -(D 21 T + D 22 T X )) (121) 

The eigenvectors of Du constitute the aforementioned diagonal basis Oi. In this frame 
D\\ — — Fx, by which token the pole part of Di is the negative of the one-loop dilatation 
operator. 

We will now consider the epsilon expansion of equation f)106p order by order in g 2 up 
to 0(e). For the rest of this section we assume the operators to be eigenvectors of -D n . 

We may take X(L, s) out of our system of equations: any set of renormalization 
factors, that renders finite the bare correlators without the X(L,s) factor, remains a 
solution on multiplication by X(L, s) because the latter is not singular in e. 

From the constant part at g° we immediately identify a Qi = t 0ii . At 0(g 2 ) the 
epsilon expansion yields simple logarithms, simple poles and a constant part. From the 
first two sets of terms and the diagonal of the third we learn 

Jn — z m — —Dun, a-ii — D Wii toa, (122) 

while the off-diagonal part of the constant term constrains B but is not sufficient to fix it 
completely; hermiticity of the two-point function on the l.h.s. of (11201) halves the number 
of independent equations. (This places constraints on D w . Similarly _D 2 i is constrained 
by the hermiticity of the l.h.s. of f)12ip .) 
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At 0(g 4 ) there is a number of conditions to solve: the double pole and the double 
logarithm in the epsilon expansion of (1 1 6 [) yield two equations implying that 

2-m = D 22il = ~ 7l 2 , (123) 

while the log(x 2 2 )/e terms give nothing new. The diagonals of the simple logarithm and 
simple pole parts lead to 

72i = z 21i = -2(D 2 i« - D 10H D 11H ). (124) 

The r.h.s. of the last equation is actually the action of a combination of Di, D 2 : 

~(D 21U ~ D 10ii D nu ) = ((D 2 - l&f)Oi). (125) 

The off-diagonal entries of the simple logarithm part depend on B and those of the simple 
pole part on B and C\. The matrix C\ cannot yet be fixed uniquely, but we now have 
enough equations to compute B. The resulting matrix equation is the off-diagonal part 
of 

B D u - D U B = -2 (D 21 - ^(D n D w + D w D n ) - \{D n D 10 - D 10 D n )) ■ (126) 

Remarkably, the last term in this expression does not contribute on the diagonal, because 
Dn is diagonal. Hence the matrix 

D 2 = -2(D 2l -±{D n ,D 10 ] - ±[D n ,D 10 \) (127) 

has 72 i on its diagonal and it determines B through (11261) . 

It was shown in [SI] that the two-loop dilatation generator acts in precisely this way: 
suppose that the dilatation operator has an expansion 

A = 1 + g 2 A 1 + g*A 2 + ... . (128) 

We want to solve the eigenvalue problem 

A(0 + g 2 BO + ...) = (1 + g 2 Y l + g 4 T 2 + ...)(0 + g 2 BO +...). (129) 

Here F\, T 2 are diagonal matrices containing the anomalous dimensions of the individual 
operators, and the lowest order re-mixing of the operators is named B. The dilatation 
operator acts on the vector of operators O as a linear map 

A 1 = D 1 0, A 2 = D 2 0. (130) 

Once again, we choose the basis for the operators to be the set of eigenvectors of D\, so 
that Ai O = D\(D = YiO. The eigenvalue problem at order g A yields 

D 2 = T 2 + {TtB - STi) (131) 

exactly like D 2 from (11271) . Note that the diagonal of B remains undetermined — it 
corresponds to trivial operator rescalings and may be put to zero. 

We have thus identified the matrix elements of the two-loop dilatation operator from 
the renormalization procedure in dimensional regularization. The next section addresses 
the construction of the dilatation operator itself. 



31 



B.2 The one-loop dilatation operator 

In the planar limit the combinatorics for the two-point functions (X y) has the following 
features: 

• At tree-level, we find a cyclic sum over, say, site 1 in X joining site % in y. All 
other lines are parallel. 

• At loop-level, the interaction is between adjacent sites. It can occur at any site in 
each part of the tree-level configuration. 

The 0(g 2 ) contribution to the correlator (Xy) originates from the M = 2 supergraphs 




Figure 3. Graphs defining the one-loop dilatation operator. 



where, of course, the underlying Feynman integral is the same in Go a , Got- It was 
called Go in [38] and is one-loop divergent. The "BPS-like" graph B Q is finite. The 
third structure F simply has a free vector line; it involves no loop-integration. The 
configurations Go a , Gq6 occur with the gauge line emanating from any of the four end- 
points; likewise, F a , F b must be joined by the opposite constellations. 

It is natural to interpret the one-loop interaction as a sum over a two-site "Hamilto- 
nian" shifting over all sites in X, which is then contracted on y much as in the tree-level 
correlator. The combinatoric factors for the Feynman graphs can be found by looking 
at the correlator 

F 1 (s 1 ,s 2 , s 3 , s 4 ) = (Tr{T a D^ZD s ^Z){l)Tr{DfZDfZT b ){2)) (132) 

at leading order in N (i.e. N 2 ), which is in a manner of speaking the one-loop interac- 
tion excised from the full correlator (Xy). We find a -2 for the "disconnected parts" 
Gob, Fb and a 1 otherwise. The disconnected diagrams can be attributed to the two- 
site interaction to their left or to their right, so that we scale by 1/2 in order to avoid 
over-counting. 

If the interaction connects sites i, % + 1 in X to j,j + 1 in y, then the other fields in 
the operators are joined by parallel free lines 

n (si , S2 ) = d:\dH n 12 = -6 SUS2 & l 2 s ^lt^ ] +■■■, (133) 

where the omitted terms contain x z or x s . The key observation is that the X(L,s) 
term in the complete correlator can only exists when all free lines have the same spin 
at both ends Coupling between sites with different spin is only possible where 

the interaction is; since we want no x± 2 with free indices the interaction can at most 
"transfer" a derivative from one of the two sites to the other. In particular, it must 
conserve the total spin. 
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Let us normalize by the inverse of the tree-level. This will simply remove all the free 
lines and scale down Fi 



Ji® — Fi(si, Sj + i, Sj, Sj + i) 



1 



n(sj, Sj) U(s j+ i, s j+ i) 



(134) 



where F\ is F\ with the over-counting corrected and the group factor N 2 5 ab stripped off. 

Whithout any derivatives, the graphs Go a , Gob, F a , Fb are absent while B = 0(e), 
whence A(0, 0, 0, 0) — > 0. When the total spin is not zero, TC[ (s) is conveniently given 

matrix: 
At spin 1 we can have 



{s u s l+1 }, {Sj,s j+1 } e {{1,0}, {0,1}} 
and our set of graphs produces 



(135) 



ttf(l) 



(136) 



At spin 2 we have the basi 



{ Sl ,s l+1 }, { Sj ,s j+1 } E {-{2,0}, {1,1}, -{0,2}} 



(137) 



and the rules for transferring derivatives are 



nf\2) 



( 2 -1 -A \ 

'2 ^ 2 * 



V 



At spin 3 the basis elements are 



i 

2 
1 
" 2 











\ 

1 

2 



{ Si ,s i+1 }, { Sj ,s j+1 } e {^{3,0}, ^{2,1}, 1{1,2}, ^{0,3}} 
while the derivatives may be transferred according to 



138) 



(139) 



/ Ai _i 



nf\3) 



-i 



^ -i 

2 x 



-1 



V4 



+ 



1 

6 / 



\ 



\ o o o \ ) 



(140) 



The pole part of these rules accurately reproduces the result of [34]: The diagonal entries 
are h(si) + h(si+\) where h(n) are the harmonic numbers, and the off-diagonal entries 
are —\jd where d counts the number of transferred derivatives. The finite part could 



17 The normalization of the basis elements reflects the fact that several derivatives at the same site 
are indistinguishable. 
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doubtlessly also be fitted: we observed that the contribution from Bq apparently always 
equals that of F a , F b which is trivial to compute. Graphs G 0a , G 6 contain only a one- 
loop integral, so that a result can be obtained in closed form. On the contrary, at the 
two-loop level this is not easy due to the complexity of the integrals. Consequently, we 
limit the scope of this work to the first few cases obtained by direct calculation. 
The one-loop dilatation operator is defined as 

i 

d 1 = Y j k 1 , (141) 

i=l 

i.e. the "Hamiltonian" runs over all sites in an operator X, mapping it to a sum of terms 
with a new distribution of the derivatives over the sites in the chain. By construction, 

((Di X) y) g o = (Xy) 9 2. (142) 

We conclude the section with two remarks: first, the definition of Ti!p in (11341) is 
necessarily asymmetric because we have normalized from the right. Correspondingly, the 
constant parts of the transfer rules are not symmetric matrices. On the other hand, the 
pole part is symmetric, because in terms of complete two-point functions the matrices 
To and I\ must be simultaneously diagonalizable. Second, it should be stressed that the 
X(L, s) terms are by far better suited to the construction of the interaction Hamiltonian 
Tii than for example the terms with no traces considered in [38J: those allow non- 
vanishing free lines between D^Z^) and Dl 2 Z(2) for unequal spins Si ^ s 2 , and the 
interaction need not conserve the total spin either. While the pole part of the one- 
loop dilatation operator is correctly obtained in this picture, we found it problematic to 
consistently subtract out disconnected parts at two loops. 

B.3 The two- loop dilatation operator 

In analogy to (11321) we try to read off the operator D 2 from the 0{g i ) contribution to 

F 2 (s 1 , s 2 , s 3 , s 4 , sb, s 6 ) = ( Tr(T a D?ZD S /Z DfZ){l) Tr{DfZ DfZD s ^ZT h ){2) ) . 

(143) 

In doing so we should remember that matrix elements of the two-loop dilatation operator 
were defined by several terms, most prominently 72 came about as a matrix element of 
the combination D 2 — Df/2, see equation (11241) and the comment after it. We fall on 
the renormalization scheme of [13J: the two- loop effective vertex has to be corrected by 
subtracting the square of the one-loop vertex. Explicitly, we take out 

\d\ = i£«f£«r (144) 

i 3 
i+Kj i 

(The derivation of the dilatation operator presented here is "asymptotic" in that it 
assumes the existence of disconnected pieces.) The first term in the last formula cor- 
responds to the situation where the two one-loop Hamiltonians do not overlap, thus all 
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terms are disconnected. If both pairs {i,i + 1}, {j,j + 1} are outside our "window" 
F 2 , they will simply cancel disconnected parts that we do not see in the excised part. 
Likewise, if only one of Ti,\ , Ttf touches the excised part, we would see an order g 2 
contribution, which we need not consider. Thus the cases of interest are (we put the left 
of F 2 at position i) 

(i) Wfw^, («) (Hi) nf} x Hfl 2: (145) 

whose relevant g 4 diagrams may be directly subtracted from the set of graphs in F 2 . 
The second term in (I144p is unfortunately not amenable to this treatment: by way of 
example we do not have a diagram that identically equals two consecutive contributions 
of G 0a . 

Our strategy thus starts by setting up an operator Jj from the g 4 graphs in F 2 
with the subtraction of disconnected parts described in the last paragraph, whereas the 
overlapping part of (_Di) 2 /2 will be dealt with later on. To avoid over-counting we have 
to rescale contributions with free lines: in complete analogy to the one-loop case we scale 
down by a factor 1/2 such graphs, that connect two matter lines but leave the right or 
left line free. Note that no re-scaling is needed when the free line is the central one; 
this situation is particular to exactly one position of the Hamiltonian. Configurations 
with two free lines can be arbitrarily shifted between the three positions within the 
Hamiltonian because the dilatation operator will involve a sum over positions. In order 
to compensate over-counting we choose to scale by 1/4 if the interaction is concentrated 
on one of the outer lines, and by 1/2 if it is on the central line. We define 

Ji = F 2 ( Si , s i+1 , s l+2 , Sj , s j+ i, s j+2 ) ] — — , (146) 

Ll{Sj, Sj) ll{Sj + i, Sj + i) ll{Sj +2 , Sj +2 ) 

with F 2 being F 2 after the appropriate modification of the set of graphs and once again 
after omission of the group factor N 3 S ab . 

The connected part in (11441) can be derived from the transfer rules for derivatives 
given in the last section. Recall that according to equation (11271) the matrix elements of 
the two-loop dilatation operator also contain the term — l/A(DuD w — D w Du), when 
the dilatation operator is made to reproduce the 0(g 2 ) remixing B. By splitting the 
one-loop transfer rules into a pole part Tif} and a constant piece H^-, we can construct 
this term as an operator in much the same way as the connected part of D\. Note that 
— 1/4 [7i^, H^j] has no disconnected part since 7^ and commute when they do 
not overlap. 

Finally, the full two-loop dilatation operator takes the form 

1 

D2 = £«?V, (147) 

i=l 

with the two-loop Hamiltonian 

nf } = ^ (148) 



2 

1 1 1 <i/o WO) ^(o WO) , «/(o) -1/(0) 1 -1,(0) -1,(0) Wo)_<i,(o)> 
^2 !' oi + n n n 0i+l + n n+i n oi + 2 0i+1 ~^ 1 ' 



4 
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The Hamiltonian Tif^ has in fact a non- vanishing 1/e 2 part, but the second order poles 
are distributed over the matrices in such a way that they drop in the sum over all 
positions. The transfer rules below and in the main text describe the 1/e part. 

In this appendix we give the transfer rules corresponding to the full Hamiltonian 
including the terms in the last line of f)148p . These come from the commutator 



-1/4 [H$ 



(0) -1/(0)1 
i > rL j l 



Note that this term is anti-hermitian, so that the transfer rules cannot be transformed 

into symmetric matrices. In the main text we omit the commutator term, since in the 

context of the Bethe ansatz it is preferable to have a hermitian Hamiltonian. In any 

case, the exact resolution of the mixing is not easy to obtain in the Bethe ansatz picture 

which projects out the descendants. Surprisingly, the formulae below do apply to the 

length 3 spin 3 mixing problem although they were derived for longer chains for which 

disconnected pieces have to be subtracted. 

The explicit bases and two-loop transfer rules up to spin 3 are: 
Spin 1 

basis: {{1,0,0}, {0,1,0}, {0,0,1}} 



V 



i 



i 

x 2 



1 - 



/ 



Spin 2 

basis: {±{2,0,0}, {1,1,0}, {1,0,1}, ±{0,2,0}, {0,1,1}, |{0, 0, 2}} 



/ 



Hf\2) 



19 

32 

21 

16 

3 

•1 

1 

4 

1 



16 

9 

4 



11 
16 



1 

"-1 



3 
2 



1 

1 

" 2 



1 

2 

29 
16 
1 
2 

67 

16 

29 

16 

1 

2 



11 
16 
9 
4 

17 
16 



•A\ 

1 

'8 
3 
•1 
1 
4 

21 
16 
19 



iy _iy i 

16 32 / 
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Spin 3 

basis: {f{3,0,0}, f{2,l,0}, |{2, 0, 1}, 2, 0}, {1, 1, 1}, ±{1, 0, 2}, 
|{0,3,0}, f {0,2,1}, ±{0,1,2}, |{0,0,3}} 
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11 
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77 
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288 
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144 
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16 
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24 




43 
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43 
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16 
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8 




91 
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n 
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11 


1 


191 
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n 
U 
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16 


1 
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29 
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16 


16 
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16 
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17 
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1 




48 


8 


4 


2 


16 


32 


3 


2 




1 


19 





19 








1031 


19 


19 




216 


144 


18 


144 


18 


144 




1 





1 





11 





191 


81 


5 




18 


4 


16 


72 


16 
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J_ 

72 

48 
J_ 
" 18 
J_ 

24 



1 

'216 
91 
144 
43 
72 
25 
288 



C Explicit Solution for Twist-Two 

The one-loop Bethe equation f|T6l) may be recast as a second-order difference equation 
for the Baxter-Q function Q s (u) 



T s (u) Q s (u) 



{u + l -) L Q s (u + i) + (u- l -) L Q s (u - i) , 



(149) 



where Q s {u) is a polynomial of degree s in the variable u, whose algebraic roots are the 
Bethe roots {uk}, 



Q s (u) = C s Yl (u - Uk) 



(150) 



k=l 



i.e. the solutions of ( ]16|) . and C s is an, for our purposes, irrelevant normalization constant. 
For twist L = 2 the excitation number s has to be even, and the Baxter equation (11491) 
is exactly solvable in terms of a hypergeometric function 



Q s (u) = 3 F 2 [-s,s + l, - -m;l,l;l] with T s (w) = 2w 2 



S "2- 



(151) 



The hypergeometric series terminates if s is an even natural number, and therefore 
generates the explicit polynomial solution of the twist-two Baxter equation"*^!. The roots 
are all real and their distribution is even, i.e. the Q s {u) in (11511) are actually polynomials 
in u 2 . Therefore the cyclicity constraint in ( TT7T) is automatically satisfied. The energy is 
found from (fTTD . (fToU]) to be 



2i 



A. 

du 



log Q s (u + 



u=0 



4 U{s+ 1) -V(l)) =4/i(s). 



(152) 



18 The details of the solution (|15ip were worked out by Virginia Dippel (unpublished) by adapting 
the method of 39J to the present case. The polynomials (| 1 5 1[) belong to the family of so-called Hahn 
polynomials [39) . 
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Here h(s) = Ylj=i Vi are the harmonic numbers, which may also be expressed through 
the logarithmic derivative of the gamma function ip(s) = d/dslogT(s). In practice, the 
roots are easily found with a root finder. E.g. with Mathematica one may define 

Hahn [s_ , u_] : = 

Expand [HypergeometricPFQ [{-s, s + 1, 1/2 - I u}, {1, 1}, 1]] 

and generate a table of all Bethe roots up to spin s, with an accuracy of k digits, 

utable[s_] := Table [Flatten [NSolve [Hahn [2 t, u] == 0, u, k]], {t, 1, s/2}] 

This is suitable, without further refinements, for finding the Bethe roots up to spin 
s ~ 70 with an accuracy, if desired, of hundreds of digits. 



D Fourier Transforms 
D.l The Gauge Theory Ansatz 

In this appendix we find the Fourier Transform of the fourth term on the r.h.s. of (172]) . 
in which is the the density u{u') is integrated against the kernel 

= % ^ l0g { l-J,2x-\u)x+\J) ) • (153) 
The definitions used in the last formula are 



« = «(«) + , M = |(l + ^-^), (!54) 

1 

— u ± - , x ± (m) = a^u 1 * 1 ) . (155) 

The branch cut of the square root is defined by the principal branch of the logarithm. 
In the following we parametrize by 

,1 . _ 1 _ . . 

u = - — iu = —iu , u = - + iu = iu , (156) 

which obey the relation 



^fWf = m ± (157) 
because both u + , u~ have positive real part. Further, let 



= a/ 1 + -V _ (158) 



and 



Kq(u) = — , Kf{u) = -^—\ 1 - — ^— ] . (159) 
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Since we are, in this paper, exclusively interested in symmetric densities, we will consider 
aw' <-> —vl symmetrized version of the kernel. Our principal equation is 



,2 



ig* / dXX 



o 



«9 U (V+H - K (u)) (k+{u') + if -(«' 
+ d u (K+(u) - K^{u)) (K+{u) + K^(u 



(1 - g 2 /2x+{u)x-{u'))(l + ^ 2 /2x+(m)x + (m')) 



2 <9« log _ g 2/ 2x -(u)x+(u'))(l + g 2 /2x-(u)x-(u')) ' ' 

To prove this, we first do the parameter integrals on the left hand side: 

d\ = ^og(uy(u) +u'y(u')), (161) 

y{u)y{u') 2g 2 

Xy(u)y(u') 

rfA-^rr = log(A)-log(l + y(u)). 

Here we rely on (j!57p to simplify. Next, we change back to the original variables u^. We 
express the roots by u, x(u) using the second relation in (11541) in the form 



1 _2^ = 2 £ M_ 1 (162) 

and finally eliminate w in favor of x(u), g 2 by the first relation in (j!54p . In a last step 
we collect all terms into one logarithm and factor the argument. As long as g is small 
this will not shift the logarithm by some multiple of 7r; one may check that the Fourier 
transform below commutes with the Taylor expansion in g. 
Next, we observe 

POO 

Kf{u) = J dte ±lut e~ t/2 Jj(V2gXt), j = 0, 1 (163) 

and hence 

/oo 
dte iut e~ W2 Jj(V2gX\t\), (164) 
-oo 

/oo 
dte iut \t\e~ ltl/2 Jj(y/2gX\t\). 
■oo 

Summing up, we have shown that 

1 ( {I- g 2 /2x+(u)x-(u')){l + g 2 /2x+(u)x+(u>)) Y 

2 u ° g \{l - g 2 /2x-{u)x+{u')){l + g 2 /2x~{u)x-{u')) 1 1 >0) 

/oo poo 
dte lut / dt'e lu,t ' \t\e-^ + ^ 2 K(V2g\t\, V2g\t'\) 
■oo J — oo 
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where 



k(t, t') 



Jo(Xt) Jo(Xt') + Ji(At) Ji(At') 



Ji(t)J (f) - J Q {t)J x {lf) 



t - V 



(166) 



We conclude that for symmetric cr(w') 



141/2 / due~ itu / du'Ktu, u') cr(u') 



(167) 



/OO P 
dt'K(V2g\t\, V2g\t'\) e^ 2 
-oo L J — 



-oo 

"OO 



-l*'l/ 2 / du'e iu,t 'a(u') 



/oo 
dt'K(V2g\t\, V2g\t'\) a(t') 
-oo 



where we have used ( I7"3l) . Now, cr(t') is an even function if cr(-u') is. We may thus reduce 
to the positive half axis, which yields the final form of the last term in (1741) . 
Similar to formula (11631) one has 



dt e ±iut e~ t/2 



V2gt 



(168) 



(V2g)i-i 



J 



u± { 1 + v 7 ! + 2£ 2 /(^) 2 



3 > 1- 



From this one can easily derive the following pretty result for the Fourier transforms 
q r {t) of the eigenvalues q r (u) of the commuting operators of the integrable magnet. The 
expression [TB] 



q r (u) 



r — 1 \ v x+(u) r_1 x~(u) 



turns intct^l 



qJt) = [ due- ltu q r (u) =4tt ( — 
J-oc \ 19 



r-2 



r-1 



-1*1/2 



(169) 



J r ^{V2gt) 

Vigt 



In particular, using this result for r = 2 we obtain the expression (1771) for the energy 
E(g) in Fourier space. As a further corrollary we find the Fourier transform of the third 



term on the r.h.s. of (172]), as stated in (1741: 



" |i|/2 r° due~ itu -d u ( 1 + — 



u 



2 ire 



.\t\ MV2g\t\) 



V2g 



(170) 
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These expressions were first obtained by Didina Serban (2005, unpublished). 
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D.2 The String Dressing Factor 

In order to include the dressing factor for the "string Bethe ansatz" [4*31 [20] , as needed in 
the discussion at the end of section 13. 3[ we replace in equation (1721) the kernel K(u, u') 
from the gauge theory Bethe ansatz by 



K s (u, u) = -d u (u - u) lo£ 



(1 - g 2 /2x + (u)x~(u'))(l - g 2 /2x-{u)x + {u')) 



(l - g 2 /2x+(u)x+(u'))(l - g 2 / 2x-(u)x-(u')) 

(171) 

which will again be needed in a u' <-> — vl symmetrized form. In view of the analysis in 
the last section, the question arises as to whether this expression can also be written as 
a one-parameter integral over pairs of the form K^(u) K^{u'). As we shall see shortly, 
this is indeed the case. 

Quite clearly we have to deal with two distinct pieces, namely the part involving 
d u u and that with d u u' . In the first case, the expression is explicitly symmetrized in u' 
whereas d u u on the whole is also even with respect to the integrals on the half axis that 
we may expect to find. We are led to look for combinations involving K^{u) + Kj{u) 
and likewise in u' . We remark that under the Fourier transform d u u <-> —td t . The 
differential operator can thus be incorporated at no expense. Surprisingly, the Kq alone 
suit our purpose: In the same fashion as before we may demonstrate 

1 ( {I- g 2 /2x+(u)x-(u'))(l - g 2 /2x-{u)x+{u')) \ , , 

O v U log -, ri r + [U <r-> —U 

2 U - g 2 /2x+{u)x + {u')){l - g 2 /2x-{u)x-{u')) ] 

POO POO 

,2 I ja. Aut I j+1 Au 1 1' 




= -g 2 dte lut / dt'e lu 1 x (172) 

J ~oo J ~oo 

2 \t\ <9| t | e-CW + l*'!)/ 2 ! dXXJ (V2gX \t\) J (V2gX \t'\) . 
Jo 

For the second piece we must try K^{u) — Kj{u) and similarly for u', because the simple 
derivative in u is odd while the extra power of v! forces antisymmetrization on the log 
factor. In a beautifully symmetric way we can realize the term as a parameter integral 
this time over antisymmetric combinations of only K\\ 

1 ( {I- g 2 /2x+(u)x-(u'))(l - g 2 /2x-{u)x+{u')) Y 

2° uU i0g l (1 - g y2x+(u)x+(u')){l - g 2 /2x-{u)x-{u')) 1 l " ' " ' 



f OO POO 



= -g 2 dte lut / dt'e lu 1 x (173) 

J —OO J —OO 

2 \t\ d W \ e -(W + l*'l)/ 2 / dX X JtiV^gX \t\) Ji(V2gX \t'\) . 
Jo 

In the right hand sides of the last two formulas the derivatives can either fall upon the ex- 
ponential or on the Bessel functions. Accordingly, we reproduce the Fourier transformed 
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gauge theory kernel K and an additional piece y/2 g K , denned as 

K(t,t') = -2 f dXX \d t J (At) Jo(At') + d# Ji(Xt) Ji(At')] (174) 
Jo L -I 

t[Ht)Ut') - Jo(t) J 2 (t')] 
(t - t')(t + #) 

Here one should first do the parametric integration in both terms separately and then 
differentiate and simplify. 

Equation (11671) is replaced by: 

/oo /*oo 
due~ ltu I du' K s {u, it') a{v!) = (175) 
-oo J — oo 

/oo 
dt' [K(V2g\t\,V2g\t'\) + V2g K(V2g\t\, V2g\t'\)] a(t') . 
-oo 
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